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We study the response of hadron masses with respect to chemical potential at p = 0. Our preliminary results of 
the pion channel show that dm/dn in the confinement phase is significantly larger than that in the deconfinement 
phase, which is consistent with the chiral restoration. 
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1. Introduction 

As suggested by QCD sum rule analysis 0, 
hadron masses may be affected by density effects. 
This may explain some results of heavy ion col- 
lision experiments such as dilepton spectra and 
J/ 5" suppression. 

It is difficult to introduce density effects in 
lattice QCD calculations due to the well-known 
"complex action" problem. Here we calculate the 
response of hadron masses to chemical potential, 
dm/d^, on dynamical configurations with /i = 0. 
Since simulations are done at fi = 0, there is no 
difficulty in obtaining dm/dii. We investigate the 
dependence of dm/d/j, with the temperature. 



where Sg is the standard plaquette action and 



Sf = ^TrlnAf(t/,M) 



(2) 



where M{U, fj) is the staggered fermion Matrix. 

The zero momentum hadron correlation func- 
tion G{t) is given by 

G(t)=^<i7(x,i)i/(0,0)t> (3) 

X 

and 

< H{x,t)H{Q,0)^ > (4) 
dUH{x, t)H{0, Oy exp{-Seff)/Z 



2. Formulation 

We use 2 flavors of staggered quarks. The ef- 
fective action to simulate Nf fermion flavors is 



Seff - Sc 



(1) 



*Talk presented by T.Takaishi 



where Z is the partition function. 

Taking a derivative of the hadronic correlator 
with respect to fi. 



d < H{x,t)H{0,0)^ > dCix.t) 

= < — > (5) 



< C{x,t)—^ > + < C{x,t) >< > . 
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where C{x, t) = H{x, t)H{0, 0)t. We calculate eq. 
(^) on dynamical configurations with /i = 0. In 
the case of /z = eq. ^ can be simplified using 
the following facts: 

(A) dSp/dfJ' corresponds to the ferniion num- 
ber operator. Thus, the average of the fermion 
number operator at ^ = is zero: < >= 0. 

(B) On each configuration the value of dSp / dfji 
is purely imaginary Thus, the value of 
< C{x, > is also purely imaginary provided 
that the operator C(a;, t) is real. This is indeed 
the case if we consider C{x,t) for mesons made 
up of degenerate quarks. 

Using the facts (A) and (B) above we derive 



d < H{x, t)H{0, 0)t > _^ dC{x,t) ^ 



for mesons consisting of degenerate quarks. 
In the spectral representation, 

Git) =J2A, cosh(m,(t - Nt/2)). 



(6) 



(7) 



Taking a derivative of eq. (0) with respect to 
fi we obtain 



dG{t) 

drrij 



Y}-pr^cosHm,{t-Nt/2)) (8) 



+ ■ 



dfi 

A{t - Nt/2) sinhK(t - Nt/2))]. 



Our procedure to obtain dm/ dfi is as follows. 
First we determine Ai and by fitting correla- 
tion function data to eq. (|^). Substituting the 
values of Ai and rrii into eq. (^) we fit the data 
of ^g^*-* to eq. (H). Then we obtain dnii/dfi and 
dAi/dfi as fitting parameters. 

3. Definition of d/dfi 

We study the two flavor case {u and d quarks). 
In this case, we have two independent chemical 
potentials, /i„ and fid- Instead, the following com- 
binations are convenient, fis = [fiu + fid)/2 and 
fiv — (flu — A*<i)/2 whith fis the usual chemical 
potential corresponding to baryon number. Then 
derivatives with respect to fis and fiy are 

d d ^ 8 d 
dfis dfiu dfid dfiu dfid 



8 8 d 8 d 
8fiv dfiu dfid dfiu dfid 

For degenerate systems of u and d quarks, 



(10) 



dfli 



dCud 

dfiu dfid 



= 



(11) 



at flu = fid = 0- In this study we analize d/dfiy 
which gives non-trivial results even with degen- 
erate quarks. In the following d/8fi stands for 
d/dfiv- 



P=5.20 m=0.025 



o dG,(x)/d^(*-l) 
□ GJx) 




Figure 1. The pion correlation function, G^(x) 
and its derivative with respect to the chemical 
potential, at /3 = 5.20. gives neg- 

ative values. To plot them in logarithmic scale, 
they are multiplied by -1. Single pole fitting re- 
sults are also shown, represented by solid lines. 



4. Preliminary results 

We present preliminary results of dm /dfi for 
Nf = 2 staggered quarks. Simulations are done 
on a lattice of size 16 x 8 x 8 x 4 at = 0.025 
with f3 = 5.20, 5.26, 5.32 and 5.34. We use the 
R-algorithm to generate configurations. The fi- 
nite temperature transition occurs at /3 w 5.28 
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Figure 2. dm/ of light-light and heavy-light Figure 3. Response of the coupling A to chemical 

pseudoscalar mesons as a function of T/Tc- potential, dlnA/dji as a function of T/Tc- 



|3| and the above (3 values are translated to 
T/Tc = 0.90, 0.97, 1.06 and 1.09 respectively. 

We measure the pion screening mass. The 
quark propagator is calculated with niq = 0.025 
(light) and 0.25 (heavy). Then we construct the 
pion correlator with light-light and light-heavy 
quarks. 

Fig. 1 shows the pion (light-light) correlation 
function GTr{x) and its derivative with respect to 
H at f3 = 5.20. We perform single pole fit for the 
data, which turned out to be sufficient for the 
pion channel. 

Fig. 2 shows 9m/ 9/i as a function of T/Tc. De- 
spite the large errors we observe a systematic ten- 
dency towards raising the derivative of m above 
T 

Fig. 3 shows the response of the coupling A, 
dhiA/dix as a function of T/T^- Both light-light 
and heavy-light mesons show similar values and 
no apreciable temperature dependence. 

5. Discussions 

Our preliminary results show remarkable char- 
acteristics of the response of meson masses to 
chemical potential. Possible interpretations for 
dm/dji of the light- light system are as follows. 



The weak response of the mass below Tc indicates 
a persistence ofthe Nambu-Goldstone boson na- 
ture at least up to T = 0.97Tc. Growth of it above 
Tc is consistent with chiral restoration since the 
meson looses the Nambu-Goldstone character. 

Calculations reported here were done on HSP ( 
NEC ) at INSAM, Hiroshima Univ. This work is 
supported by the Grant-in- Aid for Scientific Re- 
search by Monbusho, Japan (No. 11740159). 
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